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We investigate the phase transitions of holographic p-wave superconductors in 
. (4 + l)-dimensional Einstein- Yang-Mills-Gauss-Bonnet theories, in a grand canoni- 

cal ensemble. Turning on the back-reaction of the Yang-Mills field, it is found that 
the condensations of vector order parameter become harder if the Gauss-Bonnet co- 
efficient grows up or the back-reaction becomes stronger. In particular, the vector 
order parameter exhibits the features of first order and second order phase transi- 
tions, while only the second order phase transition is observed in the probe limit. We 
discuss the roles that the Gauss-Bonnet term and the back-reaction play in changing 
the order of phase transition. 
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INTRODUCTION 



> 

in 

The AdS/CFT correspondence |lj-|4j provides a novel approach to study the strongly 
\ coupling systems at finite density. Therefore, it may have some useful applications in^con- 



densed matter physics. It has been applied to study the holographic shear viscosity [5H8J, 



holographic superconductors [9|, [10| and holographic (non)fermi-liquids j 1 lM 1 31] . In this pa 



per, we will focus on the holographic p-wave superconductors 14N16 



The phenomena of superconducting can be explained by the spontaneously breaking of 



U(l) gauge symmetry [181 ] . In the p-wave superconductor, the rotational symmetry is also 
broken by a special direction of some vector field in addition. This could be achieved by the 
condensation of a charged vector field. And the holographic modeling of this picture could 
be simply realized by adding an SU(2) Yang-Mills field to an AdS black hole background [13]. 
In this case, after making some ansatz of the SU(2) field, one U(l) subgroup of SU(2) is 
considered as the electromagnetic gauge group, in addition, a gauge boson generated by an- 
other SU(2) generator is charged under this U(l) subgroup through the nonlinear coupling of 
the non-Abellian gauge fields. In this setup, the superconductor phase transition is studied, 
and conductivities show some anisotropic behavior. Further more, the holographic p-wave 



superconductor with back-reaction was investigated in Ref. [16[. The crucial point is that 
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the back-reaction will dramatically change the order of the phase transition. More specifi- 
cally, when the matter field coupling goes beyond a critical value, the former second order 
hase transition with a 1/2 mean- field theory critical exponent near the critical temperature 



., 1m, ll9N22| will be changed to a first order phase transition. 

In a previous paper [l9j . we studied the holographic Gauss-Bonnet p-wave superconduc- 
tors in the probe limit. The holographic Gauss-Bonnet superconductors are also discussed 
23-28|. In this paper we will investigate the holographic p-wave superconductors with 



in 



back-reaction in order to find out how the matter couplings and the Gauss-Bonnet coeffi- 
cient affect the phase transition of the p-wave superconductor. We find that the bigger the 
Gauss-Bonnet coupling is, the bigger the condensation value of the order parameter is, and 
the lower the critical temperature is. This reflects that the big Gauss-Bonnet coefficient 
will make the superconducting phase transition hard, which is consistent with our previous 



conclusions 19[. Besides, we also find that the stronger the matter field couples to the 



background, the harder the condensation to be formed. In addition, we find that the phase 
transition will change from second order to first order when the back-reaction is strong, 



which is similar to the discussions of [161 ] . In grand canonical ensemble, we study the free 
energy and entropy of the p-wave superconductor which also support our claim of the change 
of phase transitions. 

This paper is organized as follows: We will set up our model of the holographic super- 
conductors in Sec. [Ill and study the condensation behavior of the vector order parameter 
for different Gauss-Bonnet coefficients and different matter field couplings in Sec. IIII1 In 
Sec. IIVI we study the thermodynamics of the p-wave superconductor by exploring the free 
energy and entropy. We draw our conclusions in Sec. |VJ 



II. HOLOGRAPHIC SET UP OF P-WAVE SUPERCONDUCTORS 

We consider the Einstein-Gauss-Bonnet gravity with an SU(2) Yang-Mills field in (4 + 1)- 
dimensional asymptotically AdS space-time. The action is 

S = j ' d\^g[^(R+^ + ^{R 2 -AR^R, v +R^^ (1) 

where K5 is the five dimensional gravitational constant with 2k\ = 167rG5, and G5 a (4 + 1)- 
dimensional Newton gravitational constant, g is the Yang-Mills coupling constant and L is 
the AdS radius. The SU(2) Yang- Mills field strength is 

F; v = d,A a v - d u A; + e ahc A\A c v (2) 

where a,b,c = (1, 2, 3) are the indices of the generators of SU(2) algebra, /i, v — (t, r, x, y, z) 
are the labels of space-time with r being the radial coordinate of AdS. The A^ are the com- 
ponents of the mixed- valued gauge fields A = A^^dx^, where r a are the SU(2) generators 
with commutation relation [r a , r b ] = e abc r c . e abc is the totally antisymmetric tensor with 
£ i23 —_)_]__ The quadratic curvature term is the Gauss-Bonnet term with a the Gauss-Bonnet 
coefficient and R^ vpa = d p T ,1 ua . — ■ ■ ■ . Sbdy includes boundary terms that do not affect the 
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equations of motion, namely the Gibbons-Hawking surface term, as well as counter-terms 
required for the on-shell action to be finite. We will write S^y term explicitly in Sec. IIVI 
The Einstein field equations can be derived from the above action as 

Rp, - \aA R + ^) + f [ H i»> ~ = K l T ^ ( 3 ) 

with 

T, v = ptr(Jf^JW _ Ig^F^): (4) 

H^v = 2RR fll/ + 2R fiaTp R l/ arp — AR lip R p l/ + AR a p R p Uflcr , (5) 
H = H^. (6) 

where "tr" takes the trace over the indices of SU(2) generators. The Yang-Mills equations 
of motion are: 

V ll F afJ " u = —e abc A h F c ^ v (7) 



Following Refs. |14Hl6l| , we choose the ansatz of the gauge fields as 

A(r) = 4>(r)r 3 dt + w(r)r 1 dx. (8) 

In this ansatz we regard the U(l) symmetry generated by r 3 as the U(l) subgroup of SU(2). 
We call this U(l) subgroup as U(l)3. The gauge boson with nonzero component w(r) along 
x direction is charged under A 3 = 4>(r). According to AdS/CFT dictionary, <f)(r) is dual to 
the chemical potential in the boundary field theory while w(r) is dual to the x component 
of some charged vector operator J. The condensation of w(r) will spontaneously break the 
U(l)3 gauge symmetry and induce the phenomena of superconducting on the boundary field 
theory. 

Our metric ansatz following Ref. 16|, |17| is 

ds 2 = -N(r)a(r) 2 dt 2 + j^dr 2 + r 2 f(r)^dx 2 + r 2 f(r) 2 (dy 2 + dz 2 ). (9) 

with 



where m(r) is a function related to the mass and charge of the black hole. The reason for 
this metric ansatz is that: the back-reaction of nonzero w(r) will change the background 
of space-time. The condensation of w(r) will preserves only SO(2) symmetry of the spatial 
direction, i.e., (y, z)-direction. The horizon of the black hole is located at while the 
boundary of the bulk is at r^ y — > oo. Note that when r — > oo, 



i 



So we can define an effective radius L c of AdS space-time as 



From this relation we can see that in order to have a well-defined vacuum for the gravity 
theory, there is an upper bound for a < L 2 /4. The saturation a = L 2 /4 is called Chern- 
Simons limit. If we further consider the causality constraint of the boundary CFT, there is 
an additional constraint on the Gauss-Bonnet coefficient with — 7L 2 /36 < a < 9L 2 /100 in 



five dimensions 29-35 



The Hawking temperature of this black hole is 



T 



aN> (a . 2 0' 2 V 

t\i _ , 

r=r h vVrL 2 9 \2t\o) 



47T 



r=r h 



(13) 



where, n g = K^/g is regarded as the effective matter field coupling, and " ' " denotes the 
derivative with respect to r. The Bekenstein-Hawking entropy of the black hole is : 

s - A _ - ^ Fr 3 (u) 

S -4G 5 ~ n 2 ~ (14) 

where A denotes the area of the horizon and V — J d 3 x. 

The Einstein and Yang-Mills equations of motion with the ansatz (181191) can be explicitly 
written as 

K 2 g rf 4 w 2 (j) 2 K 2 r 3 (f)' 2 



3r 2 Na 2 - QaN 2 a 2 {rf- 2 )'(rf) 12 3r 2 a 2 - QaNa 2 (rf~ 2 )'(rf)' 2 
r(NL 2 (6f 2 + 6r 2 f 2 + K 2 g fw' 2 ) - 12/ 2 r 2 ) + 2AL 2 raN 2 f (r f)' (f^r f')' 
3L 2 r 2 f 2 - 6L 2 af 2 N(rf- 2 )'(rf)' 2 '' 

n 2 rf A w 2 (j) 2 K 2 r 3 (p /2 

a - 



(15) 



6r 2 N 2 a - 12aN 3 a(rf- 2 )'(rf)' 2 Qr 2 Na - 12aN 2 a(rf~ 2 )'(rf)' 2 

12r 3 af 2 + L 2 a (K 2 g rNw ,2 f - 3/ 2 (iVV 2 + 2N(r - a(rf- 2 )'(rf) l2 N')) + Qr 3 Nf 2 ) 
+ 6L 2 fN{r 2 - 2aN(rf- 2 )'(rf)' 2 ) ' 

(16) 

f" = h + h + f3, (17) 

? = + (-- + -)#, (18) 

r z N r a 

// ™<t> 2 .,1 J' N' a\ 

w = -TjTl ~ w - + 4 T + ^7 + - ' 19 
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where 

h = (-^rVW + ^/ 5 ^V(2/-r/0(r/)')/(3rVW 

+3ar 2 N 2 a(2r 2 Naf' 2 - 2r 2 ff'(aN' + 2Na') - f 2 (raN' + 2N(ra)')) 
+6a 2 rN 3 a(rf)' 2 (aN' + 2Na')) (20) 



h = «Jr 2 V - n 2 g afNa(2f -r/')(r/)V 2 J / (3r 4 /V 

+3ar 2 (2r 2 iVCTf 2 - 2r 2 ff'(aN' + 2Mr') - f 2 (raN' + 2iV(ra)')) 

+6a 2 rN(r/)' 2 (aAf' + 2iW) J (21) 



/ 3 = ( -L 2 r 3 «7/'(/(r«77V / + 7V(3«7 + rff / )) ~rNaf)f 3 

+raff [(L 2 a 2 N' 2 r 2 + iV (k 2 0' 2 L 2 + 2aiV'a'L 2 + 12a 2 ) r 2 + 2L 2 iV 2 a (2a + ra')) f 
+rf (-AL 2 N 2 a 2 + L 2 r 2 N' 2 a 2 + riV (K 2 r0' 2 L 2 + 2raN'a'L 2 + 2a 2 (6r - L 2 N'))) f 
-2L 2 r 2 Naf' 2 (2raN' + 3Nf (2a + ra')) - AL 2 r 3 N 2 a 2 f' 3 ] 

+AL 2 a 2 N 2 af (f + rff (-f 2 + r/7 + r 2 f 2 ) (aN' + 2iVa') ) j 
(h 2 r A Na 2 f + 2L 2 ra 2 A^ 2 a (/ + rff (aN' + 2Na') f 2 

+L 2 r 2 aNa (2Naf' 2 r 2 - 2ff (aN' + 2iVa') r 2 - f 2 (raN' + 2N (a + ra'))) ,/ 2 ^ . (22) 
There are four useful scaling symmetries in the above equations: 

(/) / ->• A/, w^A'V (23) 

(II) a^Aa, 0-+A0, (24) 

(///) r^Ar, m^A 4 m, w -> Aw, -> A0, iV -> A 2 iV (25) 

(IV) r — > Ar, m — )• A 2 m, L — >■ AL, — >• A~V, /% — > \n g , a — > \ 2 a. (26) 

We can use the symmetries ( )25|26|) to set r^ = 1 and L = 1, and use symmetries ( J23|24|) 
to set a(r — > oo) = /(r — > oo) = 1 in order to make the solution asymptotically approach 
to an AdS. 

III. TWO BRANCHES OF SOLUTIONS 

A. Analytic charged GB-AdS solution with w(r) = 

In the ansatz (jHJ) with w(r) = 0, there is an straightforward analytic black hole solution 
which is a charged generalization of GB-AdS black holes 36|, |37| . In this case, 



f(r) = a(r) = 1, <f> = A* " (27) 



6 



and 



N(r) 



r 

2a 



1 - 




A f M 1 



K 2 



9L 

9Q r 6 



(28) 



where M is related to the mass of the black hole, Q is the charge of the black hole and 
Q = 1[ir\. Using the formula ffl3|) . the temperature for this charged GB-AdS black hole is 



T 



'3nr h 



(29) 



B. Superconducting solutions with w(r) ^ 



In order to investigate the superconducting solutions with w(r) ^ 0, we have to solve 
numerically the set of equations (115 ] (115 )1 ([17)1 (ITS )) and (115 )1 . 

First of all, we should impose boundary conditions on the fields. The horizon is located 
at r = Th with N(rh) = 0. On the horizon, we should impose 0(r^) = for the U(l)3 gauge 
field to have a finite norm, and cr(r), f(r), w(r) should be finite. We can expand the fields 
in the powers of (1 — rh/r) near the horizon, they are 

' N = 0^m = rl/L 2 + m^ } (l - r h /r) H 



Near Horizon < 



a 



a 



(0) 
H 

c(0) 



a 



f =f ( H > + f ( H>a 



(1), 



rh/r) + 
rn/r) + 



(30) 



w 



to 



»8>< 

(o) 



H 



+ W 



r h /r) + 4>%\l 
(tit 



r h /rf 



H 



rh/r) + 



where all the coefficients of the expansions are constants. 

At the boundary r — > oo, the asymptotical behavior of these fields are 



Near Boundary < 



■in 



a 



f 



w 



m 



a 



(0) , 

B + m B 

(4) 



2) /r 2 + 
3 /r* + - 



f B 0) + f B 4) y + 



4 0) 

(0) 



(31) 



w 



B 



do) 



t(2) 



p, where /x and p are respec- 



(0) 



is the source 



From the AdS/CFT dictionary we know that 

tively the chemical potential and density of the charge at the boundary; w 

(2) (2) 

of the boundary operator J while w B is the expectation value of J and the nonzero w B 
will induce superconducting phase as we have mentioned. 



Near horizon (1301 ) . 



a 



(0) ,(0) 



my,a H , 



J H 



c(2) 



and w/9 can be evaluated as functions of 



H ' 



and by substituting the expansions into the equations of motion. There- 



fore, there are four independent initial values left to be specified, i.e., f B \ 
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T/Tc 




FIG. 1: The condensation values of vector operator J versus temperature for different a and 
different n g . The black, red, blue and pink curves correspond to K g = 0.0001,0.2,0.4 and 0.45, 
respectively. 

For the boundary conditions at infinity, we impose cr(r) = /(r) = 1 to have an asymptotic 
AdS boundary, this could be reached by using the scaling symmtries fl23|) and fl24l . We also 
impose = to turn off the the source of J. In the numerical calculations, we will set 
r/j = L — 1 by using the scaling symmetries (125]) and (T26l) . 

Armed with these equations of motion and boundary conditions we can numerically solve 
the set of equations by the shooting method. Because we will work in the grand canonical 
ensemble, we can fix the chemical potential value /i and then vary the four independent 
near-horizon coefficients {cr^\ fg\ 4>^h ■> w iP) until we find a solution which produces the 
desired value of \i and = fjp = 1, = 0. 

In the following, we will present our numerical results of the condensation value of the 
order parameter J. We will scan through value of a from a = —0.19 to a = 0.09 as well 
as K g from n g = 0.0001 to n g = 0.45. From the AdS/CFT dictionary we know that the 
conformal dimension of vector field in five dimension is A = 3, therefore, J 1//3 /T c is the right 
dimensionless quantity. 

Fig. [1] shows the condensation value of vector operator J for different Gauss-Bonnet 



s 



TABLE I: The condensation values of J for different a and k 9 , numbers in boldface represent first 
order phase transition 



a 

K g \ 


-0.19 


0.0001 


0.01 


0.09 


0.0001 


5.89103 


6.01942 


6.02764 


6.10352 


0.2 


6.86288 


7.35303 


7.38871 


7.74425 


0.4 


12.28375 


16.42425 


16.77913 


20.71994 


0.45 


16.17033 


24.76784 


25.65623 


35.20014 



coefficients and different matter field couplings. We can see that the condensation value 
grows if the Gauss-Bonnet coefficient grows or the matter field coupling grows. Table. [I] lists 
all the explicit condensation data of Fig. CD • 

Take the bottom-left plot in Fig. [T]as an example (i.e., the plot with a = 0.01), when we 
decrease the temperature, the condensation values for black (k 9 = 0.0001) and red (k 9 = 0.2) 
curve will emerge from zero at some critical temperature T c . When we keep on cooling 
down the system, those condensation values will continuously tend to some constant values. 
Besides, when T ~ T c the condensation will take a mean- field theory critical exponent 1/2 in 
the form of J cx (1 — T/T c ) 1//2 . These are the second order phase transitions which have been 
explored in our previous paper [19|. However, for the blue (k 9 = 0.4) and pink (k 9 = 0.45) 
lines the behavior is much different from the above two curves. We see that these curves 
will bend to the right for T > T c and there will be two condensation values for T > T c . 
However, when T < T c they will also tend to some constant values. The mean-field theory 
critical exponent 1/2 never exists for these condensations. We will argue in Sec. IIVI that 
these peculiar condensation behaviors are features of first order phase transitions and the 
real critical temperature for these phase transitions in not exactly T c . 

In the holographic p-wave superconductors, the normal charge density p n is 0^ near 

(2) 

the horizon, while the total charge density is pt = 2(f) B = 2p where the factor 2 appears 
due to the scaling behavior of <fi on the infinite boundary. The superconducting charge 
density is defined as p s = p t - p n , Es|- We plot the ratio p s / Pt in Fig. [2j For the black 
and red curves which are of second order phase transitions, the intersecting points where 
the curves meet the horizontal axis represent the critical temperatures T c /p at which the 
superconducting phase occurs. However, for the blue and pink curves which are of first 
order phase transitions (see Sec. HVj) . the intersecting points are not the critical temperature 
for the phase transitions. But the real critical temperature T C R for the first order phase 
transition is a little above the intersecting points T c . From Sec. \IV\ we can find the values 
of T^/p for a = 0.01 and varying n g s, which are listed in Table. [Til Despite of the first 
and second order phase transitions, we still see that the dimensionless critical temperature 
T^/ p decreases when k 9 grows, which reflects that the stronger the matter field couples to 
the gravity the harder the phase transition occurs. For sufficiently low temperature, all the 
curves will tend to one, which reflects the fact that the superconducting charge is dominated 
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a=0.01 

1.0 
0.8 
0.6 

pJPt 

0.4 
0.2 
0.0 

0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 

FIG. 2: The ratio of superconducting charge density and the total charge density versus the 
dimensionless temperature T j \x at a = 0.01. The black, red, blue and pink curves correspond to 
Kg = 0.0001,0.2,0.4 and 0.45, respectively. 

TABLE II: The quantity of T^/fi for different n g while fixing a = 0.01. Boldfaces represent first 
order phase transition. 



Kg 




0.0001 


0.07885 


0.2 


0.06437 


0.4 


0.02880 


0.45 


0.02063 



in the total charge. 

IV. THERMODYNAMICS 

A. Euclidean action and counter-term methods 

From the AdS / CFT correspondence, a non-extremal black hole corresponds to a thermal 
equilibrium state at the boundary. The Hawking temperature of black hole is the temper- 
ature of the boundary field. In the following we will work in the grand canonical ensemble 
with the fixed value of chemical potential /i. The partition function of the bulk theory is 

z = e -iB\g.] j (32) 

where -Te[<7*] is the Euclidean action evaluated on the on-shell value of g. Because of the 
Euclidean action, the compactified time direction has a period 1/T. The free energy now is 




Q = -T\ogZ = T I E [g*}. 



(33) 
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In the path integral, the Euclidean on-shell action should additionally include the Gibbons- 
Hawking surface term to give the correct Dirichlet variational problem and some other 
boundary counter-terms to render the action finite 38(. In the computation we introduce 
a hypersurface at large but finite r = r^y as the boundary, and then calculate the on-shell 
action by putting tm v — > oo. 

The bulk action evaluated on the on-shell values is 



rbulk 
on-shell 



V 



r 2 Na (rf)' 



f 



<*N(rf)' 
P 



2f 2 ({raN)' + ra'N) - rf (raN' + 2N (ra)') f - Ar 2 Naf 2 



(34) 



~ r bdy 



where V/T = J dtd 3 x is the volume of the (3+l)-dimensional hypersurface. The usual 
Gibbons- Hawking surface term is 



r (i) 

l GH 



h it 



d x\ 



V r 3 rrN' 



(35) 



-Tbdy 



where 7 is the induced metric on the r = r^y hypersurface, n M = \J N{r)5^' r is the outward- 
pointing normal vector to the hypersurface and K = is the trace of the extrinsic curvature 



For the GB term there is also a generalized Gibbons- Hawking term [39j, |40 



r(2) 
L GH 



■— 2 / d 4 Xy 



1 .1 - 



(36) 



aVN 



( (3raN' + 2N (a + 3m')) f - 



3r 2 f' 2 (raN' + 2N{a + ra')) f - 4r 3 Naf 



/3 



-rbdy 



where Gij is the Einstein tensor of the metric 7^ and J = ^ with 

Jij = ^{2KK ik K* + K kl K kl K i3 - 2K ik K kl K lj - K 2 K^). (37) 

No additional counter-terms for matter fields are necessary because the fields fall off 
sufficiently near the boundary 41J. 1 Besides, in our metric ansatz the scalar curvature R 
for the hypersurface is zero. So the simplest counter-term is 431446 



let — —o d £V^7— 



3 + 



V 
T4 



Nar 



3 2 + U 



r=r bdy 



Thus the total on-shell Euclidean action is 



Ie[9*\ 

Then the free energy is 

Sl = T I E \g. 



j bulk 
on-shell 



r (l) 
l GH 



T {2) + T 

GH + 1 ct- 



T(I t 



bulk , r(l) , r(2) 
on-shell ' 1 GH ' GH 



l ct 



(38) 

(39) 
(40) 



1 If one works in the canonical ensemble, the charge density p should be fixed and we should add an 
additional boundary term to the Euclidean action as A/g oc / d 4 Xy/—"ftr(n p, F^ u A a ' / ). [42 1 
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B. Phase diagrams of free energy and entropy 

In this subsection, we will discuss the numerical results of the free energy Q and the 
entropy S. In Fig. [31 we plot the free energy and entropy for some typical values of a and 

Kg. 

From the first two up-left plots in Fig. [3j we find that for a = —0.19, K g = 0.4 and a = 
0.01, Kg = 0.3, the charged GB-AdS solution (purple line) ( 1251) exists for all the temperature 
T > T c and T < T c . On the contrary, the superconducting solution (blue line) exists only 
for T < T c . But for T < T c the free energy of superconducting solution is smaller than 
that of the charged GB-AdS solution, which means that the superconducting solution is 
thermodynamically favored, compared to the charged GB-AdS solution. This represents 
that when T decreases across T c , a phase transition occurs when the bulk solution going 
from the charged GB-AdS solution to a solution with nonzero vector hair which induces 
an order parameter of the p-wave holographic superconductor. In addition, we see that at 
T = T c the purple line and blue line are not only continuous but also C 1 different iable which 
can be seen from the plots of entropy S because S = —dQ/dT, see the first two up-right plots 
in Fig. [31 The entropy for these two solutions at T = T c is continuous but not differentiable. 
So according to Ehrenfest's classification of phase transitions, these phase transitions are 
second order for a = — 0.19, K g = 0.4 and a = 0.01, K g = 0.3. This is consistent with our 
remarks in the previous section. 

From the two down-left plots of Fig. [31 i.e., plots with a = 0.01, k 9 = 0.4 and 
a = 0.01, ^ = 0.45, the behavior of free energy is dramatically different from the one 
we mentioned above, there is a characteristic "swallowtail" shape of the free energy indi- 
cating a first order phase transition. Consider the a = 0.01, K g = 0.45 case for example. 
When we decrease the temperature, entering the figure along the purple line from the right, 
we reach the temperature T « 1.2T C where new solutions appear (the blue curves represent 
the new solutions). However, for now the charged GB-AdS solution is thermodynamically 
favored because it has a lower free energy. If we keep on cooling down the system we will 
still remain in the charged GB-AdS solution until T « 1.12T C . Then for T < 1.12T C , the 
superconducting solution will have a lower free energy than the charged GB-AdS solution. 
Therefore, superconducting solution will be thermodynamically favored, compared to the 
charged GB-AdS solution in the range T < 1.12T C . So the exact critical temperature for 
this kind of phase transition is Tf- = 1.12T C . This is the reason why the quantities for 
Kg = 0.4 and k 9 = 0.45 in Table. [Ill are not the exact values which the intersecting points 
denote. This is a first order phase transition due to the non-differentiable free energy at 
T = 1.12T C . From the point view of entropy, we see that the entropy will jump from the 
purple curve to the lowest part of the blue curve at T = 1.12T C (see the down-right plot in 
Fig. [3] The entropy is not continuous which also reveals the feature of the first order phase 
transition. Therefore, from the plots in Fig. [31 we can read the real critical temperature for 
the first order phase transition as T = 1.02T C when a = 0.01, K g = 0.4 and T = 1.12T C when 
a = 0.01, Kg = 0.45, respectively. 




FIG. 3: The free energy and entropy versus temperature for different a and K g . The purple line is 
for the charged GB-AdS black hole (i.e., w(r) = 0) while the blue line is for the superconducting 
solutions (i.e., w(r) / 0) . 
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FIG. 4: Classification of phase transition versus a (—0.19 < a < 0.09) and K g (n g > 0). The order 
of phase transition depends on the coefficients a and K g . The yellow region below the dashed curve 
is of the second order phase transition while the orange region is of the first order phase transition. 

V. CONCLUSIONS 



In a previous paper 19], we studied holographic p-wave superconductors within Gauss- 
Bonnet gravity in the probe limit. In this paper, we continued this study by including 
back-reaction of Yang-Mills field. We found that both the Gauss-Bonnet coefficient and 
back- reaction will make the superconducting condensation difficult. This difficulty can be 
seen both from the growing condensation values and the decreasing critical temperatures. 
By studying the thermodynamics of the system in grand canonical ensemble, we found two 
kinds of phase transitions of the holographic p-wave superconductors. Note that in the probe 
limit, the superconducting phase transition is always second order. With back-reaction, we 
found that when fixing the Gauss-Bonnet coefficient, there was a critical value for the matter 
field couplings K 9 ( c )(see Fig. H]). If K g < K g ( c ), the phase transition is second order (yellow 
region), however, if n g > K g ( c ), the phase transition becomes first order (orange region). 
It was found that the stronger back-reaction will not only make the condensation value 
bigger but also will change the order of the phase transition. This is consistent with the 



conclusions of 16 



On the contrary, if we fixed the matter field couplings n g , 



the Gauss- 
Bonnet couplings would change the order of the phase transition just for a small range of 
K g , i.e., 0.366 < n g < 0.427. However, out of this range, the Gauss-Bonnet term would not 
change the order of the phase transition. Therefore, we may conclude that although both 
the Gauss-Bonnet coefficient and back-reaction will make the p-wave condensation hard, the 
back-reaction plays a major role in changing the order of the phase transition. 
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